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Abstract 



We study the effect of mixing two rules on the dynamics of one-dimensional cellular 
automata by large scale numerical simulations. We calculate the decay of the mag- 
netization for the Domany-Kinzel automaton (XOR/AND mixing) to its equilibrium 
value in the three phases. This requires system sizes in excess of 1 million sites. We 
also find severe finite size effects near the new critical points recently proposed on 
the basis of transfer matrix arguments. 
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ABSTRACT 

Wc study the effect of mixing two rules on the dynamics of one-dimensional cellular 
automata by large scale numerical simulations. We calculate the decay of the mag- 
netization for the Domany-Kinzel automaton (XOR/AND mixing) to its equilibrium 
value in the three phases. This requires system sizes in excess of 1 miUion sites. We 
also find severe finite size effects near the new critical points recently proposed on 
the basis of transfer matrix arguments. 



Cellular automata have been studied in a systematic way for more than ten 
years [|l| and constitute a deceptively simple family of dynamical systems-where 
space, time and the range of values of the dynamical variable have been discretized; 
yet still there are many open problems. Even better, one finds new features in well 
studied systems. 

As a concrete illustration we shall study the Domany-Kinzel automaton 0, de- 
fined as follows: 

F(afi = l|a*„i = l,a*+i = l)=p2 
P(a*+i = l\aU = 1, <i = 0) = P(a*+^ = l\aU = 0, = 1) = Pi (1) 
P(0|l,0) = l-pi; F(0|1,1) = I-P2; P(1|0,0) = 

starting from some initial configuration {af^^^f^^. If pi = 1 and J92 = one readily 
identifies the XOR rule, F(1|1,0) = F(1|0, 1) = 1, P(l|l,l) = 0, while, if = 
and P2 = 1 one has the AND rule, P(l|l, 1) = 1 all others giving zero. 

This system has been the object of many investigations since its introduction in 
ref 10-0, ^, §, §], since many models of interest may be mapped onto it, in some 
region of its parameter space. Nevertheless many open questions remain, first among 
them being the accurate determination of the transition points and a study of its 
critical properties |^. Whereas the original papers 0, Q identified one critical point 
(from a mapping to an exactly solvable model), namely at pi = l,p2 = 0.5 and 
conjectured the existence of another one (but could not resolve the issue due to the 
small system size), recent work has not only-essentially- settled this issue (in the 
positive sense) 0, ^, 0, but has advanced the conjecture of the existence of novel 
critical points within the active and the chaotic phases 0. 

In the present note we present results from large scale simulations on the decay of 
the magnetization in the various phases of the model. We also investigate what hap- 
pens near the (critical) points {pi = 0.614,^2 = 0.925), and {pi = 0.648, p2 = 0.861), 
conjectured to exist on the basis of approximate renormalization group calcula- 
tions 0. 

First of all one finds that the system sizes studied to date (even those of ref. 0]) 
are too small to resolve this sort of question. As an example we display in fig. |l| 
the magnetization, m(t) against time for system sizes up to IMsites. Especially in 
the active and chaotic phases, equilibrium is reached so fast that it is impossible to 
measure how the magnetization decays to its equilibrium value. The breakpoint (in 
the frozen phase) appears to be 128ksites. Therefore all our runs were henceforth 
done at 128ksites and larger. In fact, even this size was too small to resolve many 
situations and 4 Msites was barely satisfactory (see below). On the other hand. In 
the frozen phase we find exponential decay with characteristic time r ~ 80, fig. [l|. 

In the active and chaotic phases, we find exponential decay, with a characteristic 
time essentially equal to that in the frozen phase for the active phase; the chaotic 
phase, on the other hand, has a significantly longer decay time, as might be expected. 
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Figure 1: Magnetization vs. time in the frozen phase-(pi = 0.3,^2 = 0.7) for several 
sizes-4096,65536,131072,262144,524288,1048576 from top to bottom. 
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Figure 2: Magnetization vs. time in the active phase-(pi = 0.75, p2 = 0.87) for 
4194304 sites. 
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Figure 3: Magnetization vs. time just inside the chaotic phase-(pi = 0.76, p2 = 0.2) 
for 4194304 sites and two different initial conditions. 

Regarding what happens at the points R and Hi, we find that the values given 
for pi,p2 are almost surely misleading-the finite time effects from the transfer ma- 
trix method must be much too strong. For instance, for L = 131072 we find an 
exponential decay with definitely different decay times-cf. Fig. . 

To summarize, we have studied the time dependence of the magnetization in 
the three phases of the Domany-Kinzel cellular automaton. Our results show that 
4 Msites barely allow one to reach the asymptotic regime slowly enough to extract 
a relaxation time. Such large systems and long times are needed to resolve the 
differences (if they do exist) between the true points Hi and R. On the other hand 
it should be noted that severe critical slowing down near these critical points is a 
problem at these sizes, as preliminary simulations have indicated. Even larger sizes 
are necessary to study the properties of higher correlation functions f^-however these 
systems are definitely within reach of contemporary supercomputers and even fast 
workstations (memorywise there is no problem in running even 16 Msite systems 
on a Sparc 10 (or even an ELC) workstation and on a Cray YMP memory and 
time requirements fall within the "tiny" class). These large sizes and long times 
are necessary; since already, for much smaller sizes, we find, for instance, non- 
zero damage at pi = l,p2 = 0.38 (for L = 65536 and 1% initial damage) and for 
Pi = 0.74, p2 = 0.0, (for L = 45056 and 1% initial damage)-cf. ref. [0. 
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Figure 4: Normalized magnetization, M(t) = m{t) — m(oo) vs. time for the point 
(pi = 0.641, p2 = 0.861) (point R) and the point (pi = 0.614, p2 = 0.925) (point 
Hi); L= 131072. 
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